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SUMMARY 



The field for curved blades is represented by a vor- 
tex series with a vortex removed at the blade point. 
Further, an example of calculation of a curved blade from 
this series is given, whereby the necessary accuracy re- 
quired of the different methods in practice is shown ac- 
cording to the case considered. 



INTRODUCTION 



In the new theory of impeller wheels, airfoil lattices 
(blade grids) have an important place (see reference l,p.2l9, 
from which figs. 1 through 3 are taken); that is, for 
circularly arranged blades (centrifugal pump, fig. l) or 
for helical sections of propeller wheels (propeller pumps, 
Kaplan turbines, fig. 2). This gives directly a straight 
airfoil lattice (fig. 3) and the former can be reduced to 
one by a simple conformal transformation (£ = In Z), 

The size and shape of the blade grid obtained for 
fixed working conditions can be secured from the exten- 
sive data for single airfoils, if the following are con- 
sidered: If an airfoil is removed from out of a grid, 
the streamlines of the flow at the point are already 
Curved through the influence of the adjacent blades. If 
one puts a thin airfoil at this point, the shape . (prof ile ) 
of which coincides exactly with a streamline, obviously, 
little force is exerted on such a blade, as it certainly 
does not alter the shape of the streamline and it thus 
exerts little force on the fluid. A force occurs if the 
ai*f8i"l*diff ers from the shape of the streamline, and the 
force and distribution of it depend essentially on the 
deviation. One can only approximately assume that the 

*"Diagramme zur Bereehnung von Fltlgelreihen. " Ingenieur- 
Archiv, Bd. II, Heft 3, Sept. 1931, pp. 359-371. 
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action of two "blades is equal if their' deviation from the 
undisturbed streamlines is equal; that is, the action for 
use of this rule does not depend greatly on the initial 
form of the streamlines.* 

It can only he approximately assumed that the action 
of two blades is equal if their deviation from the undis- 
turbed stream is equal; that is, for this rule the action 
does not greatly depend on the form of the entering 
streamlines. The practice for a single blade can thus be 
used if this shape is altered fro fit the basic flow. 
Moreover, as previously seen, the basic flow is obtained 
at the point by removing a blade, but without simultane- 
ously changing the circulation distribution of the other' 
blades. 

If the blade chord is not great with respect to the 
spacing between blades, the resulting flow can be com- 
puted by replacing the blades with vortices of equal cir- 
culation. The mathematical form representing such a 
field by a vortex series is well known (reference 1, p. 
238). The formula is inconvenient for practical use, 
otherwise we should break the field into a series for 
which a vortex is missing. We- must therefore remove a 
vortex from the field, whence we get further difficulty. 
In the region of the vortex point concerned, that is, the 
region in which we .are mainly interested, the function 
approaches infinity. 

These troublesome calculations must not always be 
repeated, as the distribution of the important values in 
the region in question is given in figures 5 to 8. Be- 
cause of symmetry it is sufficient to give only one quad- 
rant. One must only notice the sign to find the quadrant. 
To facilitate their use a small supplementary drawing is 
given in the top right-hand -corner of each plot, which 
gives information as to the sign of the single quadrant. 
Figure 4 shows the arrangement of the vortices with the 
essential dimensions and symbols. 

*Exact agreement with these assumptions is not obtained. 
Besides, appreciable variation occurs if the flow in 
which the blade is placed is accelerated or decelerated. 

(Grids with pressure rise or fall, pumps or turbines.) 

(See reference 2.) Also for strongly curved flews 

(vanes with large change of direction) are variations to 
be expected. , If one is sufficiently familiar with it, 

these variations are immaterial to the satisfactory ap- 
plication of the given diagrams. 
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The vortex spacing is a, the circulation of each 
vortex T. The diagram uses tJie dimensionless ratios 
x/a, y/a, 

u* . u ft, v* = v ^ . Wn * „ Wn i $* . . i 



For the diagrams, vortex series extending farther 

r r 

than — on one side and - — on the other parallel to 

3a 3a . 

the grid are taken as "basis. The component normal to the 
grid i9 zero. In the practical case the grid may, in 
general, have any arbitrary inflow Cj , having components 
Uj and v 0 (fig. 10). For example, the "blade series 

causes a deflection of the flow so that at unit distance 
behind the grid the velocity c 2 has components u 3 and 

v 0 . (On the basis of continuity the component v Q cannot 

change.) The generalized flow is obtained from the one 
described through superposition of a parallel flow with 

velocity c Q , the components of which are u m = 

and v 0 . To obtain the desired deflection, the change in 
Ui and u 2 , the circulation about each blade must be 
considered. 

r = (uj - u a ) a (1) 

The velocity at any point then has, if the corre- 
sponding point in the diagram has the components u* = u ^ 

and v* = v #r, the components 



U l = + (u a - u 3 ) u* 



V» ■ Ti + (Uj - U a ) V 



(2) 



The pressure rise (negative = pressure fall) is 
known (reference 1,. p. 231). 
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Pa ~ Pi = - (u a 3 - u s 3 ) - p» 



(3) 



where p* is the loss in the lattice. The -pr-e-s^sixr* ex- 
erted by a blade has the components 



P x = P v o < u i - u s) 
P y = (p 3 - Pi)a 



a 



(4) 



To explain the use of the plats, we will only deter- 
mine the proper shape and placing of the blades in a 
series, when the series must produce a given effect. We 
proceed as follows: We do not at first consider the in- 
fluence of adjacent blades, which is accurate if the 
blade chord is very small compared to the grid spacing. 
(1st approximation). We then compute the change by means 
of the charts considering the influence of the neighboring 
blades. We next imagine that the neighboring blades are 
replaced by a single line vortex, as in the diagrams. 
This gives satisfactory results for most practical cases 
(2d approximation). We can immediately drop this simpli- 
fied approach and consider the chordwise distribution of 
circulation along the blade, where we can also use the 
diagrams (3d approximation). 



CALCULATION OF BLADES WITHOUT CONSIDERING- THE 
DISTURBANCE DUE TO THE ADJACENT BLADES 



The spacing a is, in general, obtained from con- 
siderations of construction. The inflow velocity is given 
by its components v 0 and u a . The blade series will 

increase the velocity m 1 to u 3 (turbine blades). In- 
stead of the velocity change the force in the grid direc- 
tion (tangential force) P x can be used. From the gov- 
erning relation (u 3 - Uj) can be easily computed. (The 

pressure drop - p 2 ) or the force perpendicular to, the 

grid is required; thus the pressure loss p' must be es- 
timated next on the basis of corrected airfoil data.) 
The circulation r is given by equation (l) from (u a - ilj) 
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In aeronautical literature are given values of 6 a 
(lift coefficient) and e w (drag coefficient,) depending 
on the angle of attack a for wing profiles. Besides 

c w 

the so-called glide angle, the ratio = c plays a 



a 



further part in the loss p T , or the machine efficiency, 
The circulation T is related to the lift coefficient 
c a of a "blade in the lattice by 



(5) 



where t is the "blade chord in the direction of flow. 
The hypothesis for this formula is that the velocity is 
not essentially changed in the region from the leading 
edge to the trailing edge of the "blade. If this condi- 
tion is not satisfied, the "blade chord must "be corrected, 
whereupon we return to figures 11 to 13. We will now 
select a suitable airfoil from those tested, for example, 

a plate with small camber, ca-m-ber = 0.05 (reference 3, 

chord 

p. 407.* Figure 9 gives test values for this airfoil. 

In case cavitation or other considerations do not 
dictate a very low value of c a , we choose c a to give 

a satisfactory small value of the glide angle (low drag). 
In this case** for the airfoil chosen: c a = 0.7 and 

a = 3°. 

With this lift coefficient, and the magnitudes of 
Uj, and u 3 given by the problem, the necessary chord is 

fixed by equations (l) and (5). 



t_ k 2(u! - u a ) 
a »a c o 



* The values given there must be corrected to infinite 
aspect ratio. (Formulas given in Hiltte, p. 402.) 

'*Ior pressure drop through the grid, one can probably 
pick. a higher value of c a for this case. (K. 
Christian!, see footnote p. 2 of this report.) 
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For example, from the problem we require v 0 = 2u x - . aad 
u 2 = 3uj thus, 



V ' u-i + Up V 

x - = u 3 - u a = 2u a u m = 2 = 3u a = 4- 



> (7) 



°o 



= J*n 3 + v o 3 = u l ./2 s . + 2 3 = 2.82 u a = 1.41 - 
and thereby 

; - oj : 8 a . g3 - 2 - 0s (8 > 

If — is very small, we can consider the flow in the 
a 

region of the blades as uniform (1st approximation), and 

this flow has two components u m = u * + u 3 and v Q . 

2 

The deflection velocity for the blades considered at the 
zero point of the diagram leave only the superimposed 
parallel flow. We must therefore arrange the blades with 
the proper angle of attack relative to this flow. In our 
example the superimposed flow has the components u m = 2uj 

and v Q = 2u a , thus its direction forms the angle 

P = tan" 1 = 45° 'with the grid direction. The blade 
u m 

'chord has an angle of attack of 3° with this direction; 
it thus makes an angle with the grid direction of 
8 s p - a = 42°. The arrangement of the blades is shown 
in figure 10 according to the calculations. But it is 
intentionally emphasized that for the example the assump- 
tion that the -blade must be small compared to the grid 
spacing, i» not satisfactorily fulfilled, and for that 
reason the arrangement shown in figure 10 is not satis- 
factory. 



NACA Technical Memorandum No. 1022 



7 



INFLUENCE OP REPRESENTING- ADJACENT BLADES BY VORTICES 
Calculations "based on Deflection Velocities 



If the blade is not small compared to. the grid spac- 
ing, the tail thus projects noticeably into the disturbed 
region of the adjacent blade. The flow direction there 
is different from that in the middle and we must adjust 
the blade curvature for this changed relation. Por this 
calculation we can, in general-, replace the adjacent air- 
foils by separate vortices (2d approximation). Since the 
actual blade occupi es . considerabl e s"pace,we must fix a 
place in the blade to put the replacing vortex. It has 
been shown that the least error is made if the vortex is 
placed at the so-called" center of pressure, that is, the 
point at which the resultant force cuts the blade chord. 
It lies, in general, at a distance sj. between 0.25t and 

0.5t from the leading edge. The moment coefficient c m 

about Jrhe leading edge is usually given in the reports of 
airfoil research . (fig, 9). Whence one can calculate 

£..|S. (9) 
* °a 

lor our example c a = 0.7 e m = 0.29, thus 

*. - °' 29 
t ~ 0.70 



= 0.41. In figure 10 the x-axis is drawn through 



the profile center of pressure 0.41 t from the leading edge, 
The deflection velocities can be easily found for each 
point of the blade by use of diagrams 6 to 8. If the de- 
flection is not very large.it is enough to find the rela- 
tion to correct the blade at the leading and trailing 
edges, especially the latter. As the deflection becomes 
less, the position of the leading and trailing edges, com- 
puted from the deflection, approaches those for the basic 
undisturbed stream. Por v a considerable displacement of 
the edges, it can be obtained by repeated ..calculations, 
but, in general , -this wi\ll not be greatly changed from 
the assumed position. 

In our example the edges for the undisturbed flow 
lie on a straight line making an angle 6 * 42°. . The 
leading edge is removed from the zero point a distance 
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Si = 0,41t = 0.89a (t = 2.03a" from equation (7)), and the 
trailing edge a distance s 2 = 0.69t = 1.20a. 

The angular deviation of the flow follows from the 
deflection velocities u and v (charts 6 and 7) from 
the adjacent vortices. According to figure 11 

AP = tan AP = v cos B - u sin g (l0) 
c 0 + u cos 3 + v sin P 

It should he noticed that u and v have different 
signs in the different quadrants, hut P must always he 
less than 180°. For the leading edge of our example u 
and v -is- negative. In many cases the deflection com- 
ponent u cos P + v sin P«c 0 . In which case the cal- 
culations can he simplified if the deflection is not cal- 
culated for the "blade itself, hut' rather for the point 
on the undisturhed streamline c 0 through the zero point 

making an angle P with the x-axis (fig. 10). This de- 
viation component v cos P - u sin P is identical vrith 
the component w n given in chart 8, and we get 

AP ~ ± ^S- (11) 



c 



o 



The positive sign applies to the entrance distance y > 0 
and the negative sign to the exit y < 0. For our example 
we read from the chart (P = 45°, s a = 0.83a) w n = -0.14 £ 

and for the trailing edge (P = 45°, s 3 = 1.19a) w n = 0.22 
Since c 0 = 1.41 — (see equat ion (7)) , we get for the cor- 

Si 

rection angle, 

o 



for the leading edge AP « 7-| 
for the trailing edge AP ** -7 




(12) 



We must consider that the hlade does not give an appreci- 
ahle change in the velocity c 0 . For a given wing the 
circulation is proportional to the velocity. To maintain 
constant circulation we must therefore reduce this 
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velocity in inverse ratio to the "blade speed and similarly 
ITeTea s e" "it^f o r -'l'ower*>«'peed8 i .^-JLhe^e,lj3city_at the lead- 
ing edge in the example, terms in equat i on" ( 9") , 

cb = 0.74. .c 0 
and at the trailing edge. 

c A = 1.26 c Q 



If we now assume that the velocity varies approximately 

linearly from the zero point to the edges, then the mean 

c 0 + Cj; c 0 + Cj^ 
velocities for these distances are ^ or 



and from the original equations, these velocity changes 
require changes in length of the separate "blade elements 
at the inlet or exit edges of the distance Sj or s a , 

respectively, from the zero point. Where 



s' 3 = s 3 



2 Cp 
c o + C E 



3 c 0 
c o + c A 



= 1.15 



= 0.885 



Sj = 0.96a 
s 3 = 1.05a 



(13) 



With these corrections we get the changes in the airfoil 
shown in figure 11, the proportional change in the tail 
point "being clearly shown. Practically it is almost 
meaningless compared to the change in angle of the lead- 
ing and trailing edges. It is noticed "besides that for 
this change only the change in total chord of the airfoil 
has effect. If, therefore, as in the above case, the end 
points are displaced equally, there is scarcely any dif- 
ference. To "be sure this applies only when the streamline 
along which the profile lies is in some degree straight,. 
Tor strongly curved streamlines (large velocity change 
in the lattice), the displacement of the profile itself p 
plays an •important part. 

Calculation "by Conformal Transformation 

If, as in our example, the simple formula (ll) does 
not satisfactorily calculate the necessary "bend in the 
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"blade, it is then simple and accurate to obtain the curve 
by means of chart 5, in the following way (figs. 12 and 
13). 

If we draw the streamlines of the undisturbed flow 
superimposed on the deflected flow (fig. 5) to an equal 
scale, then streamlines of the undisturbed flow are 
straight lines, for our example lines at 45° to the x-axis. 

p 

Since c 0 = 1.41 — , thus choosing the distance between 

a 

these streamlines and numbering them so that streamline 
"aero" goes through the zero point and streamline "one" 

is displaced a distance — a — figure 12 shows for the 

1 . 41 

example the streamlines for which the stream function in- 
creases by 0.1. The superposition of both flows, the un- 
disturbed c Q (fig. 12, dotted) and the deflection flow 

(fig, 12, thin solid) can thus be composed by known meth- 
ods to get the resultant flow (fig. 12, solid). Similarly 
the potential lines can be obtained. Since, in general, 
only a few points need be taken, it is simple to determine 
them by vector addition. 

At the top of figure 13 the blade is shown for the 
undisturbed flow with its streamline and potential-line 
grid. If because of the adjacent blades the grid is dis- 
torted (fig. 13, bottom), the corresponding point of the 
blade, that is, leading or trailing edge lies at the 
point having the same potential and stream function (con- 
formal transformation). 



INFLUENCE OF FINITE ADJACENT BLADE DIMENSIONS 



For these calculations we have calculated the influ- 
ence of the deflection resulting from replacing the ad- 
jacent blades by a single vortex placed at the center of 
pressure of the blade. In most cases this approximation 
is completely satisfactory. Now if the blade chord is 
large compared to the spacing a cos P, it should be 
noted that the space distribution of circulation in the 
adjacent blades must be Considered. The vortex distribu- 
tion along a blade chord is at least approximately known 
from lift and pressure distribution. Figure 14 shows the 
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JJ.ft di strJLbution for the identical circulation distribu- 
tion f^- j about our«profile. It the profile is distorted 

by conformal transformation, each part is itse'lf deformed. 
The circulation about an element of the profile remains 
unchanged. Therefore, the circulation for each unit 
length must diminish proportionally to the stretching and 
conversely, as the profile is deformed. Because of this 
displacement of the circulation, the center of pressure 
is itself somewhat shifted. The computations can then 
be carried out by shading the entire profile; the change, 
however, ie relatively unimportant as shown by figure 17 
(thin broken and solid lines). On the basis of the as- 
sumed circulation distribution given in figure 14, the 
new distribution and center of pressure given by the dis- 
tortion is shown in figure 15. We proceed on the basis 
that the concentrated vortex at the center of pressure 
can be replaced by an appropriate vortex at each element 
on the blade. The influence for the vortex series at 
each element can then be separated and summed. Instead 
of working with finite elements, we can change to the 

&r 

differential ds , for which the circulation is ds 

o s 

and find the total deviation by integration along the 
chord. The latter method will be used in the example. 

There is, in general, only a very small correction 
to the calculations carried out for the concentrated 
vortex. It is thus almost completely satisfactory to 
compute only the correction for the leading and trailing 
edges of the profile obtained from equations (9) and (13). 
To illustrate the method for obtaining the calculations 
at the trailing edge, we choose about 4 points for the 
blade (a, b, c, d, fig. 15) and</etermine from charts 6 
and 7 the effect that a vortex series, of unit circula- 
tion at thiB point, has for the time being on the veloc- 
ity at the trailing edge. If the point has coordinates 
i and ti with respect to the aero point of the airfoil 
in question, and the trailing edge, the coordinates x, y, 
we thus find from our chart the disturbance due the vor- 
tex series at the point x ~ t y - r | ( fig< 16 ) t tf e thus 

a a 

obtain the two components u and v of the disturbance 
velocity at the trailing edge of the -blade, (fig. 15, 
curve E u for the u component and E y for the v 
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component) and multiplying the value of the distribution 
curve "by these influence lines, we thus get two curves 
which give the deviation due the distributed vortex at 

the blade point (fig. 15 E u |£ and E v |j £). The 

total deflection at the trailing edge of the blade dis- 
tributed vortex is obtained by. planimetering the curves, 

} u = 0.34 £■ 

v = 0.14 £ 

The concentrated vortex at 
from charts 6 and 7, 

u« = 

v 1 = 0.15 

EL 



the center of pressure gives 

0.43 £ 
a 



For the last values equation (9) gives for the angular 
def lecti on, 

v 1 cos P - u' sin P 1 o 

*»' = c n + u' cos Mt' sinp -- 6 ^ (14) 



as in figure 13. On the basis of the distributed vortex 
we get 

A f> = U 8in P - V C ° S P - - 4*° (15) 

c Q + u cos P + v sin 0 

tfe must therefore, for the vortex distributed along the 
blade chord, consider an increase in the exit angle of 

- A0' = - 41° + 6i° = lf° 



In general this small angular difference will be given 
with sufficient accuraoy by the following approximate 
equation 
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Vo Aot Xy, ..s-jeI.) cos, J - (u - u 1 ) Bin P . . 

c 0 + u' cos P + sin P 

The dl aturtance component in the small direction c 0 for 
the concentrated vortex is 



u' cos 0 + v 1 sin P = 0.29 c Q 
and for the distributed vortex 

u cos 0 + v sin 0 = 0.24 c Q 



(17) 



We must therefore displace outward the trailing edge in 
the direction c 0 from that given by equation (13) by a 

di stance 

is A = 0.02a (18) 

The corresponding calculation for the leading edge 
gives a change in angle of 

&P - AP' = -f° (19) 
and a displacement from the zero point of a distance 

Asjj * 0.025a (20) 

The profile with these changes is shown in figure 17. 
For comparison the profile from the assumed concentrated 
vortex is shown in thin lines and the profile for changed 
center of pressure due to distribution is shown in thin 
dotted lines. Thus the effect of distribution is seen to 
be fairly unimportant. 

The example was intentionally chosen so the various 
influences had discernable effect for the various practi- 
cable applications illustrated by the diagrams. In general, 
it is soon obvious which method must be used to obtain the 
profile with sufficient accuracy. Excessive accuracy for 
many cases has little sense as the profile properties for 
strongly curved streams and especially for large pressure 
drop or increase through grid are themselves not accurate. 
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(See footnote, p. 2, K. Christian!.) In many cases a 
simple angular correction from equation (11 ) is suffi- 
ciently accurate. 

The curved blade can .naturally be obtained if the 
effect of a blade grid of given shape is estimated. The 
deviation velocity is obtained, from an estimated cir- 
culation about the blade, from the shape and angle of 
attack of the blades in an undisturbed flow. By repeat- 
ing these calculations for the circulation, agreement 
with the desired value can be quickly obtained. 



Translation by J. W. McBride, 
Massachusetts Institute of 
Technology. 
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Figure 2.- Cylindrical inter- 
section through 
propeller wheel. 



Figure 1.- Airfoil arranged in 
circle. 
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Figure 4.- Explanation of symbols. 

w- velocity at point xy; 
u,v,w n components of same. 
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Figure 3.- Airfoils ar- 
ranged in 
straight line. 




Figure 6.- Velocity component u*=u^ parallel 
to lattice direction. 
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Figure 7.- Velocity component v*=v^ normal to 
lattice direction. 1 




Figure 8.- Velocity component w*n~ w nf^ normal to radius. 
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Figs. 9,10,11,12 




Figure 9.- Lift, drag, 
and moment 
of an arc shaped, cam- 
bered plate of infinite 
aspect ratio. 



Figure 10.- Form and setting of air- 
foils in first approxima- 
tion taking no account of the effect 
of neighboring airfoils). 



Figure 11.- Effect of dis- 
turbance veloc- 
ities on the profile shape 
(second approximation). Light 
line shows profile of first 
approximation. 



Figure 12.- Superposition 
of disturbance 
"flow (thin lines) on undis- 
turbed flow (dotted lines) 

fives the disturbed flow 
heavy lines ) . 





Figure 13.- Conformal transfor- 
mation of profile in 
the undisturbed flow (top) on the 
profile in the disturbed flow 
(bottom) (second approximation) . 




figure 14.- Approximate distribution of the 
circulation over the airfoil in 
undisturbed flow. 
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Figure 16.- Position of vortex at trail- 
ing edge of airfoil. 




Figure 15.- Distribution of circulation 

over disturbed airfoil and 
its effect on the trailing edge. 



Figure 17.- Profile shape with vortex 
distribution over the air- 
foil taken into account (third approx.). 
Thin line shows profile of second approx. 
(according to Fig. 13) for comparison. 
Dotted line shows position of profile af- 
ter displacement due to displacement of 
center of gravity. 



